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O ' Abstract 

^ ^ In this work, we present a supersymmetric extension of the quantum spherical model, both in components 
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3 , and also in the superspace formalisms. We find the solution for short/long range interactions through the 

Qs ' imaginary time formalism path integral approach. The existence of critical points (classical and quantum) 

is analyzed and the corresponding critical dimensions are determined. 
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I. INTRODUCTION 

It is a common fact that the methods and ideas of a given area of theoretical physics may be 
useful in another completely distinct context. Besides the technical utility that this sharing of 
knowledge may offer, it provides different views of the problem yielding profound implications. 
A classical example of this sort of situation is the junction of ideas coming from Gell-Mann-Low 
quantum electrodynamics and Kadanoff's block-spin ultimately leading to the renormalization 
group, useful both in quantum field theory and statistical mechanics [1,]. Proceeding along these 
lines of thinking, in this paper we consider the supersymmetric extension of a quantum version of 
a traditional model in statistical mechanics, namely, the spherical model [2]. Supersymmetry, by 
nontrivially combining internal and space-time symmetries in a way impossible in the traditional 
^ ^ie algebraic approach, has produced new relevant insights in the context of quantum field theory 
3(1 and even in condensed matter physics J-lSl]. In the nonrelativistic context, supersymmetry was 



introduced by Nicolai 



m 



and applied to the study of spin systems. 
In quantum mechanics, supersymmetry requires the existence of supercharges Q and Q, such 
that {Q,Q} = Ti, where 7i is the Hamiltonian of the system, and furthermore satisfy [Q,T-L] = 
[Qj^.] = 0; these supercharges realize the transmutation of bosonic states into fermionic ones and 
vice- versa so that the ground state is left invariant. Thus, a supersymmetric theory is characterized 
by a bosonic ground state lO)^ with energy equal to zero, i.e., annihilated by the supercharges and 
by the Hamiltonian 

Q\o)B = Q\o)B = n\o)B = o. (1) 

On the other hand, if the supersymmetry is broken, then, at least one of the supercharges does 
not annihilate the ground state any longer. Instead, there is a pair of degenerated ground states 
bosonic and fermionic, \0)b and \0)f, with energy Eq, such that 

|0),. = ^0|0>B and \0)s^J^Q\0)p. (2) 

In contrast with phase transitions at finite temperature, driven by thermal fluctuations and 
denoted as classical phase transition, a quantum phase transition occurs at zero temperature and 
is caused by quantum fluctuations connected with Heisenberg's uncertainty relations. Thus, at zero 
temperature, a phase transition may occur in a broken supersymmetric situation when the ground 
state energy Eq is non- vanishing; it is characterized by some non-thermal coupling parameter g, 
that assumes the value gc at critical point. So, near the critical point, 

Eor^\g-gcr, (3) 



where z and v are the dynamical and correlation length critical exponents, respectively. This 
suggests that in a situation where super symmetry holds, the system will not exhibit a critical point 
because of the vanishing of the ground state energy independently of the value of g. On the contrary, 
it may display a critical behavior if supersymmetry is broken. It should be noticed that, essentially 
because the differences in the bosonic and fermionic distribution functions, supersymmetry is always 



broken at finite temperature |9l4lll]. 

The classical spherical model was initially proposed by Berlin and Kac [2^], as a simplified con- 
tinuous version of the Ising model. Since it is exactly soluble, it has been used to study the critical 
behavior in a variety of situations 
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19|]. Typically, the classical solution exhibits an anomaly 



at low temperature, concerning to the third law of the thermodynamics. It was suggested that 
such pathological behavior could be corrected by introducing quantum fluctuations, which were 
not taken into account in the classical case [201]. The extension to the quantum domains has been 
considered by various authors and has raised much attention in the context of quantum phase tran- 



sitions 



m 



24 
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23l |. In particular, we mention studies in some quantum versions of the spherical model 
29|, and also including som e ing redients of the statistical mechanics, such as the influence of 



random fields |30|] , spin glasses 
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35l |. frustration [36|], competing interactions [37[, and quantum 



Lifshitz point [38]. In this work we extend these studies by considering a supersymmetric version of 
the spherical model with especial attention to the existence of critical points and the determination 
of critical dimensions. The supersymmetric constraints are implemented through delta functions in 
the partition function and by the saddle point evaluation we obtained the conditions to which the 
action is stationary. The critical behavior analysis follows from the study of these equations in some 
simple cases, according to the values of the saddle point parameters. For a certain situation, where 
the bosonic and fermionic degrees of freedom are decoupled, we found that the critical behavior of 
the supersymmetrical model reduces to that of the quantum spherical model. On the other hand, 
when the bosonic and fermionic degrees of freedom are coupled, corresponding to another choice 
of the saddle point parameters, the model exhibits a distinct critical behavior. We also discussed 
some issues in connection with the supersymmetry breaking. 

Our work is organized as follows. In the section |II1 we review the bosonic quantum spherical 
model and evaluate the partition function through the imaginary time formalism. In the section 
mil we present the supersymmetric extension, both in components and also in the superspace 
formalisms. The partition function of the supersymmetric model is evaluated in the section ITVl and 
its critical behavior is discussed in section |Vl A summary and additional remarks are presented in 
the Conclusions. One appendix to study the canonical quantization of the model is included. 



II. THE QUANTUM SPHERICAL MODEL 

In this section, we discuss the quantum spherical model subject to the strict spherical constraint, 
corresponding to a canonical ensemble. The quantum version of the model, by implementing the 
mean spherical constraint, was studied by Vojta [27]. The classical Hamiltonian of the spherical 
model is given by 

'Hc = \Y. Jr,r'S,Sr' + /^ ^^ ^r, (4) 

r,r' r 

where r and r' are lattice vectors, {S'r} is a set of spins variables that can assume continuous values, 
— cxD < Sr < oo, in a D-dimensional hypercubic lattice; Jr,r' is the interaction energy that depends 
only on the distance between the sites r and r', Jj-j.' = J(|r — r'|), and h is the external field. We 
assume that the S'r variables are subject to the spherical constraint 



Y,S'r=N, (5) 



where A^ is the total number of lattice sites. On the other hand, the mean spherical constraint 
is defined as J2r i'^r) = -^j where (• • • ) designates a thermal average. Of course, in the thermo- 
dynamic limit, A^ —7- cxD, these two constraints yield the same results, as we shall see shortly. A 
discussion about the properties of spherical and mean spherical models in the classical context is 
found in jl2| . 

In order to construct the quantum model we need to add to the Hamiltonian a kinetic term 
involving the conjugate momentum variable to S'r, denoted by P^, and then promote such variables 
to operators satisfying the usual canonical commutation relations 

[Sr, Sr'] = 0, [Pr, Pr'] = and [Sr, Pr'] = i^r.r'- (6) 

Including then a kinetic term quadratic in the momenta, the quantum Hamiltonian becomes 

^ = ^5 E ^r + ^ E JrySrS,, + /l J] Sr. (7) 

r r,r' r 

The coupling constant g measures the relevance of the quantum fluctuations and the limit 5 — )• 
corresponds to the classical regime. To construct the supersymmetric version it is convenient to use 
the Lagrangian formulation. Moreover, it is also useful for the evaluation of the quantum partition 
function through the imaginary time formalism, as we will explicitly show in the next subsection. 
After a Legendre transformation, we obtain the Lagrangian 

£ = — 2^Sr--2^ Jr,r'SrSr' - /l 2^ Sr, (8) 

r r,r' r 



where the dot means a derivative with respect to the time. For simphcity, in the remaining of this 
paper, we will take the external field equal to zero, h = 0. It must be clear that we can consider 
the external field dependence without difficulties. 

A. Imaginary Time Formalism 



Let us evaluate the partition function through the path integral imaginary time formalism 
approach [39]. Thus, we need to pass to the Euclidean imaginary time r = it, with r G [0, /3] and 
(3 the inverse of the temperature. Furthermore, the bosonic variables are required to satisfy the 
periodic boundary condition 5r(0) = 5'r(/3), which gives rise to the discrete spectrum Matsubara 
frequencies ujn = 2mr/f3, with n G Z. The partition function is given by 



■^ r 

where the Euclidean Lagrangian is 



(9) 






dr 






(10) 



The integration measure T>Sr symbolically stands for functional integration over the spins variable 
of all sites of the lattice, i.e., VSr = Hr ■^'^r- Employing the functional integral representation for 
the delta function 



J^52-iV) = fvXe-^od-^{^rS?-N) 



(11) 



we can write 



VSrVX exp < 



dr 






(12) 



This is an appropriate representation, since the integral on Sr becomes Gaussian and can be 
integrated out. Before continuing, however, it is convenient to introduce the Fourier transformation 
of Sr 



^j;e^.-r5,(r) 



5r(r) 



With this, the Euclidean action acquires a simple form and the partition function becomes 



(13) 



Z = / VSqVX exp 



/3 



dr 



E^. 



2g dr^ 



+ ^ J(q) + A ) 5_n - AiV 



(14) 



Here, we have identified J(q) as the Fourier transformation of the interaction energy J^ 

J(q) = ^J(|h|)e^'i-^, with h = r-r'. 
h 

After the Gaussian integration and using the identity det A = e^''^'^ , we get 



where we defined 



5, 



eff 



z 



^ 1 

^^^+2iV^^ 



V\e 



'NS. 



eff 



El- 






(15) 



(16) 



(17) 



In the thermodynamic hmit, A^ — )• oo, we may use the saddle point functional method to evaluate 
the partition function. The saddle point condition corresponds to the effective action above be 
stationary: 



6X{t) 



0. 



(18) 



We will suppose that the saddle point A is time independent, and defined as A = /i. By means of 
the identity (5Trln^ = TrA^^SA, we can obtain 



1 



^ oo 



1 



2iV/3 



q n=— oo 2g "'" 2 '^*^-' "•" ^ 

The sum on integers can be calculated thanks to the identity 



0. 



E 



vr 
^-- — 2 = - coth(7ry), y > 0. 
n"^ + y y 



The final result is 



1 



-T. 



-^coth^'^"'' 

2a;n 



0, 



(19) 



(20) 



(21) 



where, (-i^^ ^^(/x + J(q)/2). This result is exactly that obtained with the mean spherical con- 
straint 23]. From this expression, we can determine the critical properties of the model at finite 
temperature as well as at T = 0. The analysis of the critical behavior can be done by considering 
the system near the critical point, with /i — )• and the interaction parameterized as J(q) ~ q^ , 
q = |q|, for small momenta. The parameter x determines the short or long range character of the 
interaction. Typically, for short range interactions we have x = 2. Despite being described quan- 
tum mechanically, for any finite temperature the system shows a classical critical behavior similar 
that obtained from the classical version of the model. Essentially, this is because the thermal 
fluctuations generally dominate quantum fluctuations at macroscopic scales. At zero temperature. 



however, there is a quantum phase transition characterized by new critical exponents 



27[. 



III. SUPERSYMMETRIC EXTENSION 

In this section, we shall construct the super symmetric version of the spherical model in terms 
of the component fields formulation as well as in the superspace. The fundamental ingredient in 
this extension is the introduction of additional fermionic degrees of freedom at each lattice site to 
balance the bosonic ones. Namely, at each site of the lattice, besides the bosonic variable S^, we 
will associate the fermionic counterparts V'r and ipr (V'r = i^l), that in the quantum case satisfy 
the anti-commutation relations 

{'lpr,1pr'} =0, {i'r,1pr'} = and {V'r, V^r'} = f^r.r'- (22) 

A. Components Formulation 

The natural form for the super symmetric Lagrangian is 

^Susy = ^ XI "^r - 2 X] Jr,r'SrSr' + ^^ AA " ^ ^r.r'V'rV'r' • (23) 

r r,r' ^ r r,r' 

Of course, the requirement of invariance under supersymmetry transformations will imply a relation 
between the interactions Jr.r' and Uj-y. We are supposing that the interaction f/r.r' also depends 
only on the distance between the sites r and r'. In addition, the constraint on the bosonic variable 
will imply others constraints involving the fermionic variables, which will be discussed shortly. 
Now, it is easy to verify that the set of transformations 

6eSr = ■^rE) '^eV'r = ;;= Syt — N^ UrySr'e and S^'lpr = (24) 

Vs J.' 

and 

5eSr = eipr, ^elpr = and (JeV'r = —p Sre — ^^ [/r,r''S'r'e (25) 

leave the Lagrangian (1230 invariant up to surface terms, i.e., up to a total derivative, provided that 

Y,Ur,sUs,r'=Jr,r'. (26) 

s 

Notice that the Fourier transformation of this relation furnishes [?7(q)]^ = J(q). The parameters 
of transformations e and e are anticommuting infinitesimal quantities. Equations (j24p and (j25p are 
the supersymmetry transformations, that relate the bosonic degrees of freedom to the fermionic 
ones. 



The next step is to investigate the consistency of the constraints. More precisely, the imphca- 
tions of the supersymmetry transformations on the spherical constraint 



J2^r=N. (27) 



It can be verified that, under the transformations (j24p and (j25p . we are led to the additional 
constraints 

y tprSr = 0, y^ IprSr = and ^^ V'rV'r = " /^ UrySrSr'- (28) 

r r r r,r' 

In the last relation, we have discarded a surface term. In summary, to have a consistent super- 
symmetric formulation of the spherical model we need now four constraints. These constraints 
introduce effectively an interaction between the bosonic and fermionic variables. They can be 
implemented via four Lagrange multipliers, being two of bosonic character and two of fermionic 
character. 

B. Superspace Formulation 

A very elegant and concise way to formulate supersymmetry is through the notion of the super- 
space. Moreover, this formulation has the advantage of making the underlying theory manifestly 
supersymmetric. The price is that it is necessary to introduce auxiliary (not dynamical) bosonic 
degrees of freedom. The superspace consists of an extension of the ordinary space, that in the 
quantum mechanic case corresponds to the time coordinate only, in order to accommodate the 
anticommuting coordinates 6 and 6, satisfying 9^ = 9^ = and {6, 6} = 0. To fix our notation, let 
us define some useful operations with these Grassmannian variables: 

and 

f d9= f d9 = 0, I d99= I d99 = l and I'd9d9 99 = l. (30) 

We then associate at each site of the lattice a superfield $r(i) ^) d) that is a function of the super- 
space coordinates and can be expanded in powers of 9 and 9 in the following way 

$r(t, 0, 9) = Sr + Hr + i^rO + 99F^. (31) 



The bosonic and fermionic variables are components of the superfield and Fj. is the auxihary field. 
Next, let us introduce the supercharges, 

that are the generators of translations in the superspace, 

e ^ e + e, t^t-iOe and 9^6 + 1, t^t + ieO. (33) 

Also, it is easy to verify that the superfield transformations, 

5e^r = Qe^r and 6e^r = eQ^r, (34) 

correspond to the following component transformations 

5eSr = 'ipr^, (^eV'r = —iSr€ + Fj.€, (J^'i/'r = and 5f,Fj. = iipj-e (35) 

and 

6iSr = eipr, ^eV'r = 0, JfV'r = iSr^ + Fy-l and S^Fy. = —ieipr- (36) 



These transformations may be compared with (|24p and (j25p by rescaling the variables Sr — >■ 
g~^'^Sr, ipr — )■ g~'^''^'il^r and -ipY -^ g~^ 'ipr and also the parameters e — ;■ g~^''^e and e — )• g~^''^e. 
Furthermore, the supercharges satisfy the algebra 

{Q,Q} = 0, {Q,Q} = and {Q,Q} = -2i^. (37) 

The last anticommutation relation is proportional to the generator of time translations that must 
be identified with the Hamiltonian. 

We need to construct an action in the superspace such that the corresponding Lagrangian 
reproduces (j23p after integration over 9 and 9, and further take in to account the constraints (j27p 
and ()28p . Initially, let us consider the kinetic term. For this purpose it is convenient to introduce 
the covariant derivatives 



that satisfy 



{D,Q} = {D,Q} = {D,Q} = {D,Q} = and {D,D} = 2i-^. (39) 



These anticommutation relations imply that the covariant derivative of a superfield has the same 
property under supersymmetry transformations as the superfield itself. So, any action involving 
only superfields as well as covariant derivatives of superfileds is manifestly supersymmetric. 
Now, observe that 



D^r = -V'r + [iSr - FA 9 - 111)^99 (40) 

and 

i)$r = -V^r - (iSr + F^\ 9 + i4)^99. (41) 

Then, the kinetic term can be obtained according to 



2 



r r r r 



Notice that we do not have terms involving time derivative of the variable F^- , which means that it 
is an auxiliary variable as we said above, i.e., it does not possess dynamics and can be eliminated 
via its equation of motion. The interaction terms can also be constructed in a simple way 

-^ Uj, j,i^r^ri =y^Urr'SrFy,/ -S^Urr''4^r1pr'- (43) 

2 ^-^ ' ee ^-^ ' ^-^ 

r,r' r,r' r,r' 

The interaction energy Jj-j./ does not appear explicitly in the action, but only after the elimination 
of the auxiliary field F,.. The action in the superspace is then the sum of the kinetic and the 
interaction parts 

S= f dtd9d9 I - Y^ D'^rD'^r + — E ^r.r'^r^r' ) • (44) 

The equation of motion for the auxiliary variable Fy. is given by 

Fr = -Y,Ur,r'Sr', (45) 

r' 

that can be used to eliminate Fj.. The resulting action of this process is, up to a rescaling of the 
fields Sr —J- g^^''^Sr, ipr -^ 9~^'^'ipr and ipr — > 9~^''^'ipr as mentioned before, just that of components 
formulation. 

The spherical constraint will be imposed on the superfield according to 

^$r$r = A^, (46) 
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which in components yields the following relations: 

"^S^ = N, Y^ ^rSr = 0, Yl ^prSr = and Y V'rV'r = Yl ^^^'^- (^^) 

r r r r r 

After using the equation (j45|) . we get (|27p and (|28p . We can implement the constraint through 
a delta function (inside the partition function) that in the integral representation requires an 
integration over a superfield 

E = j + 9c + ce + eex, (48) 

where 7 and A play the role of the usual Lagrange multipliers and C and ^ of anticommuting 
Lagrange multipliers. To sum up, we have constructed a consistent supersymmetric extension of 
the quantum spherical model in components as well as in the superspace formulation. The next 
step is to evaluate the partition function, which will be discussed in the sequel. 

IV. IMAGINARY TIME FORMALISM 

We will employ the imaginary time formalism to the evaluation of the partition function, as it 
has been discussed in the purely bosonic model. Because of the anticommuting character of the 
Grassmannian quantities, the fermionic variables are required to satisfy the antiperiodic boundary 
conditions. In general, the fields must satisfy 

Sr{0) = Sr{/3), Vr(0) = -V'r(/3) and i;,{0) = -M(3). (49) 

The reflex of the antiperiodic conditions is the arising of the discrete spectrum fermionic frequencies 

io^ = (2n + l)7r//3, with n E Z, in contrast with bosonic frequencies uj^ = 2n7r//3. 

A. Components Evaluation 

Let us consider the partition function in components, generalizing the procedure of the pure 
bosonic model discussed before. The partition function is then given by the functional integration 
over all fields presents in the Lagrangian taking into account the constraints (|27|) and (|28|) : 

Z= fvSrVi;rVi;rd{ Y Sr-N)6{ Y i'rSr)S{ Y i'rSr)S{ Y i'vA + Y Ur,r'SrSr')e- ^0 '^^^^ ^ 
r r r r r,r' 

(50) 
where Ce is the Euclidean version of ([25|) . 

(51) 



Csusy = ^^[-Q^) + 2 2^ Jr,r'SrSr' + ^Yl ^^^ + Yl ^r^r'AA 
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As before, it is convenient to use the integral representation for the deha functions: 



5{Y,ij,S,)= VCe 



■Io^drj:^iirSr<; 



(53) 



5(^^r5r)= VCe-fo'i-T.jASr 



(54) 



and 



SiY.i'rA + Y.UrySrSr') = fvje-^odri{j:.^ry^r+j:,yU,yS.S,,) 



(55) 



In this way, the integration over Sr and (V'riV'r) becomes Gaussian and can be performed. Firstly, 
let us concentrate on the Sr integration by considering the effective partition function 



^eff ^ / -^"^^ ^^P 



dT 






dT 



)+ — y ^ JryJrOr 



+ -^ X/ "^r + 7 E ^r,r' 5'rS'r' + y^ ifrSr 



(56) 



where we defined the real bosonic field (pr = V'rC + CV'r- Introducing the Fourier transformation of 
the fields involved and performing the Gaussian integration, we find 

1 f" 



^e'}/=exp|-^TrJ]lnOqlexp|-i^''dr J^^qCOqiCV'qi, 



(57) 



with 



Or 



1 ^2 1 

+ -J(q) + A + 7t/(q). 



(58) 



Now, we move on to the fermionic integrals by defining another effective partition function whose 
action is written in terms of Fourier transformations 



^e// = / ^V'q^^q exp <^ 

which gives 

^eff = exp 



/3 



dT 



^'^'^(7^^+^^''^+^^^'^^'^'^''^'' 



(59) 



Tr y In 



J_d_ 
VgdT 



f/(q)-7-^COa'C 



1 
2CC, 



(60) 
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So, by putting all together, we get 

Z= [vW-fVCVCe-^^-ff, (61) 

with 

S.„ ^ -^^rA + J^TYX;i„(-i^ + ij(q) + A + W(c)) 

->E>"(-;;^^-c'(i)-'-^CO-c). (62) 

Proceeding as before, we use the saddle point method, which in this case implies into the following 
four conditions 

^^eff ^ SSeff ^ 5Seff ^ 5Seff ^ 

5X{t) 57(t) 5C(t) dCir) ^ ' 

We shall look now for solutions of this equations with all parameters (A, 7, C, C) time independent. 
The conditions of extremum with respect to the fermionic parameters Q and Q, can be immediately 
satisfied if (" = C = 0. So, there are still two remaining conditions. Let us define the corresponding 
saddle points as A = ^ and 7 = 0. The condition of extremum with respect to the A can be worked 
out in the same way that led us to the result (l2T]) : 



^i:^->.(9^]-. (-) 



but, with the bosonic frequency slightly different, namely, ((^q )^ = 2(7(^ + aC/(q) H — 2)' Fin^-Hy) 
the last condition furnishes 

where the fermionic frequency is (Wq )^ = 2(7(^ + aC/(q) H — ^). To obtain this result, we have 



used the identity 

00 
E f9„. ^ ^^2 ^ „2 = ^ ta^h (^) , (66) 



(2n + 1)2 + y2 - 2y ^^""^ ( 2 ) ' 



in the calculation of the trace of the fermionic part, which involves a sum on the fermionic Mat- 
subara frequencies w^ = (2n + l)7r//3, n G Z. The critical behavior of the model can be determined 
by analyzing these two saddle point conditions near the critical point, what will be done in the 
next section. 

One last noteworthy remark before closing this section is that the procedure above described 
can be also generalized to the superspace, i.e., directly in terms of the superfields. We will not 
consider this approach in this work. 

13 



V. CRITICAL BEHAVIOR 

In order to investigate the critical behavior we need to consider the system near the critical point, 
when /i — )• and a — )• 0, and the integrals are dominated by the small momenta contributions. 
Thus, as mentioned before, we can parameterize the interactions for small values of |q| = q as 
J(q) ~ q^ and U{q) ~ q^ , respecting the supersymmetry requiring [C/(q)]^ = J(q). Next, we will 
analyze the critical behavior at finite temperature as well as at zero temperature, since the system 
exhibits different behaviors in these two cases. 

A. Finite Temperature 

As already mentioned, the supersymmetry is incompatible with temperature, i.e., at finite 
temperature the supersymmetry is broken. In this situation, the thermal fluctuations are present 
and are responsible to drive the phase transition. At the critical point, the equations (j64p and ()65p 
become 



l-ly /- cothf ^V^^^^ Uo (67) 



and 



N^2./^Ji^) '''^ V 2 ; N^2./^J{^) ^"^^ [ 2 ) ' 

respectively. These integrals converge if D > x, what defines the lower critical dimension. 

We may determine the critical behavior of the system by subtracting the expression (j64p near 
the critical point from (167p . Technically, near the critical point we can expand the hyperbolic 
functions coth and tanh for small values of the argument, according to what has been discussed 
above. 

From now, we are going to investigate some simple cases that exhibit interesting critical be- 
haviors: 1. a = with finite ii near the critical point; 2. n = with finite a near the critical 
point. 
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1. a = and finite /i 

By subtracting the equation (j64|) (with a = and small //) from (|67p . and evaluating the sum 
over the momenta, we obtain the following behavior 

H^ {D < 2x) 
tg ~ < //In/x {D = 2x) , (69) 

jj {D> 2x) 

rn 

where tg = {g — gc)/gc, as in [27]. Equivalently, of course, we could consider the distance from the 
critical point given in terms of the temperature, t-p = (T — Tc)/Tc. Clearly, the critical dimension 
is given by Dc = 2x. This is exactly the same result of the purely bosonic spherical model of the 
section |TI1 We interpret this behavior as the decoupling of the bosonic and fermionic degrees of 
freedom when the saddle point value of the parameter 7 is zero (a = 0). In this situation, we end 
up just with the bosonic model subject to the spherical constraint. 

2. iJ, — and finite a 

For /x = and small a, after subtracting the equation ()64p from ()67p and then evaluating the 
remain sum over the momenta, we find 

' 2{D~x) 

a ^ 



tg ~ <; 



(^<f) 

alna (^ = f ) • (70) 

(Z?>f) 



a 



Here, we no longer have the decoupling of the bosonic and fermionic degrees of freedom, since the 
correspondent constraint effectively introduces an interaction mediated by the parameter 7, whose 
saddle point value a is now different from zero. Consequently, the model exhibits a distinct critical 
behavior and a critical dimension Dc = 32;/2. 

It should be mentioned here, despite the break of the supersymmetry because of the temper- 
ature, we could try to investigate the case when the frequencies become equal, namely, when 
/i = a^/2 (in the case of zero temperature this would correspond the supersymmetric case). The 
behavior arising from this choice is exactly that of (j70p . This can be understood due to the 

X 

dominance of the term q2 over q^ at small momenta, whenever a is non- vanishing. 
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B. Zero Temperature 

In the case of zero temperature, (3 -^ oo, the hyperbohc functions coth and tanh become 
identicahy one, and the equations (j6l|) and (j65|) are given by 



^-^E^ = » Pi) 



q q 



and 

q q q 1 

The integrals converge ii D > x/2, what defines the lower critical dimension in the quantum case. 
The procedure for determining the critical behavior is the same of the finite temperature. Note 
that at the critical point {fj, = a = 0) the equation ([72]) is identically satisfied. 

1. Supersymmetric case: 11 = ^ 

The supersymmetric situation is characterized by the equality between bosonic and fermionic 
frequencies, u^ = u^ = Wq, which can be achieved by choosing ^u = ^, as we said. The con- 
sequence of this choice is that even outside (near) of the critical point we get a = and also 
/i = 0. This result indicates the absence of the critical behavior, in agreement with the argument 
presented in the Introduction. It has a simple interpretation: the bosonic and fermionic quantum 
fluctuations necessary to drive the quantum phase transition are canceled between themselves in 
the supersymmetric situation. 

2. Broken supersymmetry case: a = and finite /z 

In this situation the frequencies are no longer equal and the supersymmetry is broken. By 
proceeding as before, we may determine the quantum critical behavior 





( 2D-X 
H 2x 


(D<f) 


tg ~ < 


/i In fi 


(D=f) 






(Z)>f) 



(73) 



where now tg = {g — g^)/g^, with g^ the critical value of g at zero temperature. The critical 
dimension is D^ = 3x/2. As already argued in the finite temperature case, this result also is the 
same of that obtained in the purely bosonic spherical model. 
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3. Broken supersymmetry case: /i = and finite a 

In this last case, where the supersymmetry remains broken, we get the fohowing critical behavior 



tg ~ < 



^U — X 

a X (D < X) 

a In a {D = x) , (74) 

a {D > x) 



which, in virtue of the coupling between the bosonic and fermionic degrees of freedom, it is different 
from ()73p . and reveals the critical dimension D^. = x. 



VI. CONCLUSIONS 

In this work, we constructed a consistent supersymmetric extension of the quantum spherical 
model, by considering the components and the superspace formulations. Afterwards, we calculated 
the partition function through the imaginary time formalism, yielding to the saddle point condi- 
tions. From these, we studied the critical behavior of the model in some simple cases for finite and 
also for zero temperature and determined the critical dimensions. In general, the model exhibits a 
critical behavior whenever the supersymmetry is broken, either by thermal effects (classical critical 
behavior) or because of the inequality between the bosonic and fermionic frequencies (quantum 
critical behavior). 

Specifically, we verified that when ^u is finite and a = 0, the critical behavior (T ^ and 
T = 0) is the same of the pure bosonic quantum spherical model, namely, Dc = 2x for T ^ and 
Dc = 3x/2 for r = 0, which agrees with the results of [27]. In the other case, when a is finite 
and /i = 0, we obtained a reduction of the critical dimension values: Dc = 3x/2 for T 7^ and 
Dc = X for r = 0. When fi = a'^/2, we have a supersymmetric situation at zero temperature. In 
this case, the model does not exhibit a quantum phase transition because of the vanishing of the 
ground state energy and consequently the absence of the quantum fiuctuations. 

As final remarks, we stress that some further issues can still be analyzed, for example, the 
calculation of thermodynamic quantities. Moreover, the supersymmetric model is an interesting 
laboratory to investigate some typical ingredients of the statistical models, such as competing and 
disorder interactions, including finite-size effects, the presence of random fields, the existence of 
quantum Lifshitz points, and their possible influences on the critical properties. Lastly, it could be 
worth to explore the connection with the supersymmetric field theoretic non-linear sigma model 



4CH42||. according to the Stanley's ideas |43|. 
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Appendix A: Canonical Quantization 

In this appendix, we consider the supersymmetric mean spherical model, with the constraints 
imposed as thermal averages, corresponding to a grand canonical ensemble. From the beginning, in 
order to simplify the analysis we will invoke the law of large numbers in the thermodynamic limit, 
what avoids the presence of orthogonality constraints in the Lagrangian because of the vanishing 
of the mixed two point functions 

^Y.i'rSr^{4^rSr)=0 and ^ ^ V^r^r ^ (V^r^r) = 0. (Al) 

r r 

In the path integral approach, these conditions correspond to the saddle point choice (, = (^ = 0. 
Without this simplification the diagonalization procedure will involve a transformation of coordi- 
nates such that Ax'^ + Bxy + Cy'^ — t- Ax"^ + By^ + C, where x and y generically represent bosonic 
and fermionic coordinates. 

Next, for our proposals it is convenient to write the Lagrangian (j23p in a more symmetric way, 
already taking into account the Lagrange multipliers 

C = Y^^^ ~2^ Jr,r'SrSr' + —p^ V'rV'r " 2 ^ ^'''''' t^^' ^^'^ 



r,r' " r r,r' 



(^^S^ - n) - a(^^Y.i'^r,A]+Y.^r,r'SrSr') - h^Sr-^{Av + fiA), (A2) 



r 



where h and (ry, rj) are bosonic and fermionic external fields. By introducing the momenta 

P, = ^, Ur = ^ and nr = ^, (A3) 

we can easily determine the Hamiltonian by means of a Legendre transformation 

r r,r' r,r' r 

+ »U^[i^r,^r] + Y.^rySrSr')+h^Sr + ^ii^rr] + fii'r)- (A4) 
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From now on, we interpret the variables as quantum operators. After writing the Hamiltonian in 
the Fourier space and then performing the standard Bogohubov transformations, we can find the 
diagonal form in terms of creation and annihilation operators 

with the frequencies (w^)^ = 2g{fi + aU{q) + ^) and {u^f = 2g{^ + aU{q) + ^) and the 
bosonic and fermionic operators satisfying 

[aq,«q'] = 0, [a^,aq,] = 0, [aq,a^,] = (5q,q/ (A6) 

and 

{Cq,Cq/}=0, {cJ,,C^,}=0, {Cq, C^, } = 5q,q/ . (A7) 

The partition function can be straightforwardly evaluated according to Z = Tre~'^ and the 
connection with the thermodynamics is through the free energy 



f = -W^''^ = -^I^ 



B 



In cosh — — — In sinh 



q 
Finally, the constraints are implemented as 






^ = and -^ = 0. (A9) 

The resulting conditions coincide with equations ()64p and (j65p . when h = and fj = rj = 0. 
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